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Abstract 

Low-scale gauge-mediated supersymmetry(SUSY)-breaking (GMSB) models with 
gravitino mass rn^n < 16 eV are attractive, since there are no flavor and cosmologi- 
cal problems. In this paper, we thoroughly study the collider signal in the case that 
the next-to-lightest SUSY particle is the bino or slepton and investigate the discov- 
ery potential of the LHC. Our result is applicable to a wider class of GMSB models 
other than the minimal GMSB models and we pay particular attention to realistic 
experimental setups. We also apply our analysis to the minimal GMSB models with 
a metastable SUSY-breaking vacuum and we show, by requiring sufficient stability 
of the SUSY-breaking vacuum, these models can be tested at an early stage of the 
LHC. 



1 Introduction 



Supersymmetric (SUSY) standard model (SSM) is the most promising candidate for the 
physics beyond the standard model (SM). Among proposed models of SUSY-breaking 
mediation mechanisms, the gauge-mediated SUSY-breaking (GMSB) model [1] is very 
attractive, since there are no SUSY flavor problems. In addition, if the gravitino (G) is 
lighter than 16 eV, there are no cosmological gravitino problems [2] . Such a light gravitino 
plays an important role in collider physics, e.g. at the LHC Once SUSY particles are 
produced at the LHC, they successively decay into lighter particles and finally into the 
next-to-lightest SUSY particle (NLSP). The NLSP then decays into a gravitino, which 
is the lightest SUSY particle (LSP). In a typical GMSB model, the NLSP is the bino- 
like neutralino or righted-handed slepton, since their gauge interactions are weak. Their 
dominant decay modes are Xi ~ > 7 + G and + G, respectively. Thus the LHC signal 
is multi-photon+missing energy, or multi-lepton+missing energy. For both signals, there 
are tiny SM backgrounds. Thus, the low-scale GMSB scenario is easier to be tested at 
the LHC compared to other SUSY models. 

There are many studies for the discovery potential of the LHC for low-scale GMSB 
models at the LHC [31 HJ El El Ej. Although these studies are basically based on the 
minimal GMSB models, general GMSB models have wider parameter space. In this 
paper, we investigate the LHC discovery potential for more generic parameter space than 
the one of the minimal GMSB models. We perform a comprehensive analysis using all 
known discovery modes and pay particular attention to realistic experimental setups. 

From a theoretical viewpoint, such a light-gravitino scenario tends to suffer from a 
serious problem: instability of the SUSY-breaking vacuum. There is a strong upper 
bound on the SUSY-particle masses for the minimal GMSB models in which the SUSY- 
breaking vacuum is metastable [8]. We apply our analysis of the discovery potential to 
such GMSB models and find that these models can be tested at an early stage of the 
LHC. Especially, if the number of messengers is one, it is possible to test at a very early 
stage: 1 fb -1 at y/s = 7 TeV. 
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2 Model Setup 



In this section, we illustrate the setup of the models we consider in this paper. 

We first consider a simple GMSB model, where a SUSY-breaking field S couples to 
N 5 pairs of messenger chiral superfields, \I/ and iff, which transform as 5 and 5 under 
the gauge group SU(5)gut- The simplest form of the coupling of the messengers and the 
SUSY-breaking field is 

W = J^iyiS + M^i, (1) 
i=i 

where Mj is the messenger mass and %)i is a Yukawa coupling constant and assumed to 
be less than 0(1). By assumption, the SUSY-breaking chiral field S develops an F-term 
vacuum expectation value (VEV) (S) = 6 2 F, which is related to the gravitino mass as 
|F| = y / 3m 3 / 2 Mp L , assuming that the SUSY breaking is dominated by F. The condition 
that Mj > yiF should be satisfied, otherwise the messenger scalars are tachyonic. 

In the minimal GMSB models, the MSSM gaugino masses are generated from loop 
diagrams of the messengers. At the one-loop level, the gaugino masses are given by 

M a = ^A 9 (l + 0(x*)), (2) 

where we have defined 

A^Ef. (3) 

1=1 * 

and Xi = ViF/Mf. Here, a(= 1,2,3) labels the SM gauge group U(l), SU(2) and SU(3), 
and we use the normalization cci = 5q;em/(3 cos 2 0\y). The scalar masses arise at the 
two-loop level, and are given by 

< = A £0 2 ^)( 1+ « ( 4 ) 

a 

where 

1=1 



Mi 



(5) 
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and C a (i) are the Casimir invariants for the visible particles 0, (C\(i) = 3Yj 2 /5). □ Each Xi 
is bounded as x« < 1 for the messengers not to become tachyonic, and then the corrections 
of 0(xf) are typically small and we ignore these corrections in the following analysis. We 
see that ~ M a = 0(1) TeV is realized for A s ~ A g = 0(100) TeV. To realize 
77i3/2 < 16 eV, the mass of the messenger must be Mj = 0(100) TeV. 

While conventional studies on LHC physics of GMSB models have been performed 
with the above setup of minimal GMSB models, general GMSB models exhibit different 
patterns of mass generation. To be more precise, in the minimal GMSB models, the two 
scales A g and A s have a certain relation. However, in many examples of GMSB models 
the relations for A g and A s are modified. For example, if the mass matrix of the messenger 
is complicated, the relation between A g and A s is also complicated [9]. Another example 
is the strongly interacting messenger. If the messenger has large anomalous dimension, 
the relation between A g and A s is deformed. Therefore, we treat the A g and A s as 
independent free parameters in order that the analysis can be applied to a wider class of 
low-scale GMSB models. 

In the most general setup of a GMSB model [10J, the pattern of mass generation is 
further more general than the one in our setup. However, essential part of our analysis 
can be also applied to this general case. We will come back to this discussion in Sec. 13.41 

Masses of MSSM Particles 

Since the above expressions Eqs. @ and (@J for the SUSY-particle masses are given at the 
messenger scale, one should solve the MSSM renormalization group equation to obtain the 
on-shell masses and mixing matrices. To calculate the on-shell masses and other physical 
parameters, we use the program IS A JET 7.72 |TT] slightly modified by the authors. We 
adopt 200 TeV as the messenger mass. 

The mass spectrum and therefore the low-energy phenomenology depend on the pa- 
rameter tan (3, which is the ratio between the up-type and down-type Higgs VEVs. In 

^"We assume that the gaugino and scalar masses are parametrized by one parameter, A g and A s , 
respectively. In general case, this universality is not maintained since the GUT is broken at the messenger 
scale and there are two types of the messengers, lepton and quark-type messengers. However, once we 
assume the both masses and yukawa couplings of the lepton and quark-type messengers are identical at 
the GUT scale, such universality is maintained. 



4 



most part of our analysis, we use the values tan /3 — 10 and 40 to represent the low-tan f3 
and high-tan cases. The effect of varying tan is discussed in Sec. HI 

In Fig. [TJ we illustrate the masses of SUSY particles. For the case with tan/3 = 10, we 
extend the range of A s to negative values, where we define "A^" to be sign(A s )|A s | 2 . This 
is because generically we may well have negative scalar masses at the messenger mass 
scale, e.g. with the D-term contribution. Since the parameter region is not excluded 
either theoretically or experimentally, we include the region for completeness. For A g , 
the lower value is experimentally excluded as we illustrate in later figures. The green 
region in each figure is the region in which the correct electroweak symmetry breaking 
cannot be achieved. In the left figures (a) and (c), the masses of the squark, gluino, and 
lightest chargino are illustrated. As we see below, these particles are important for the 
SUSY production at the LHC and their masses often determine the cross sections and 
thus discovery reaches. The masses of stau and the lightest neutralino are drawn in the 
right figures (b) and (d). These are the lightest SUSY particles in the models and we also 
draw a boundary curve where the NLSP changes. In the right side of this boundary, the 
neutralino is the NLSP and in the left side the stau is the NLSP. Thus, this boundary 
corresponds to a change in searching channel, from diphoton + missing to multi-lepton 
+ missing. 

3 LHC Signature and Discovery Potential 
3.1 LHC Signature 

In low-scale GMSB models, the gravitino plays an important role at the LHC. The pro- 
duced SUSY particles successively decay into the NLSP and finally, the NLSP decays into 
the gravitino. The decay length of the NLSP is roughly given as 



Therefore, if the gravitino is light the decay of the NLSP occurs promptly. In the present 
setup, the NLSP is the bino-like neutralino or the righted-handed slepton. Their main 
decay modes are Xi ~ > 7 + G an d $■ — > £ + G, respectively. Thus the LHC signal is 
multi-photon+missing energy or multi-lepton+missing energy. 




(6) 
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Figure 1: Contour plot of the masses of the MSSM particles. In the left figures (a) and 
(c), the squark, gluino and wino masses are shown. In the right figures (b) and (d), the 
lighter stau mass and the lightest neutralino mass are shown. Figures (a) and (b) shows 
masses for the case with tan/3 = 10 and (c) and (d) for tan/3 = 40. The green region 
in each figure shows that a correct electroweak symmetry breaking cannot be achieved 
there. In figures (b) and (d), the blue line shows the points where the lighter stau and 
the lightest neutralino have the same mass. 



6 



3.2 Analysis 



In the following analysis, we set the gravitino mass 7713/2 = 16 eV. We have used the 
programs ISA JET 7.72 [11 J to generate the MSSM mass spectrum and decay table and 
Herwig 6.510 p2] to generate SUSY events at the LHC For detector simulation, we have 
used AcerDet 1.0 [32] slightly modified by the authors. 



3.2.1 Detection Efficiency and Misidentification 

In the fast simulation, the detection efficiencies of a photon or a lepton which passes a 
certain isolation criteria are 100%. However in multi-lepton or photon signals, the reduc- 
tion of the number of signal events originated from misidentifications of these particles is 
important. In our simulation, we include the fake rates of the leptons, photons and jets. 
The fake rates we have used in the analysis are shown in Table [TJ 



Table 1: 
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ake rates. 



3.2.2 SM Background 

In the present model, high energy leptons or photons accompany the SUSY events at 
the LHC. Thus, there are few backgrounds from the QCD events. The main background 
comes from ti and gauge boson production events. To estimate the background, we have 
used the programs MC@NLO 3.42 [14J (for ti,WW,WZ and ZZ), Alpgen 2.13pj (for 
Wj, Zj and W/Z + bb/tt), MadGraph 4.1.44 [16] (for tt/W/Z + 7/77), and Pythia 6.4 
CU (for 77)- 

In contrast to the SUSY signal events, we assume that photon and lepton detection 
efficiencies are 100 % for the SM backgrounds (see Tabled]). The misidentified jets from the 
SM backgrounds can be significant for multi-lepton and multi-photon modes. Therefore, 
for the fake rate from jets to leptons and photons, the SM backgrounds are treated in 
the same way as the signal events. Therefore, our estimates for the event number of SM 
backgrounds are conservative. 



7 



3.2.3 Event Cuts and Optimization 

To illustrate the discovery potential of the LHC, we calculate the optimized significance 
of each model point against the SM backgrounds. To incorporate both the statistical and 
systematic uncertainties of the backgrounds, we adopt the method used in Ref. [5]. The 
systematic uncertainties of the backgrounds are taken to be 50% for the QCD multi-jet 
background and 20% for others. 

Given the number of signal events N s , background events N b with the uncertainty SN b , 
the significance is given by calculating the convolution of the Poisson distribution with 
some "posterior" distribution function. As the posterior distribution, we take the gamma 
distribution as suggested in Ref. [18]. The resulting significance Z B is given by [18] 



Z B = v / 2err 1 (l-2p B ), (7) 

with 

B(N S + N b ,l + Nl/Ml 5N 2 J{N b + 5N*)) 
PB B(N S + N b ,l + N^/6N^) ' w 

where erf -1 is the inverse error function and 

B(a,b,x)= dtt a ' 1 (l-t) b - 1 (9) 
Jo 

is the incomplete beta function and B(a,b) = B(a,b,l) is the usual beta function. If 
we take the limit 5N b — > 0, the Eq. ([8]) reduces to the probability in the usual Poisson 
distribution. When N b — 0.1 and N s < 10, Zb — N s . In the case of smaller background 
N b < 0.1, we conservatively take the N s as the significance. 

For each model point, significances, as defined above, are calculated with some sets of 
kinematical cuts and a set of cuts which gives the largest significance is selected. After 
this optimization, the condition Zb > 5 is used as a criterion for discovery. 

3.2.4 Search Modes and Kinematical Cuts 

As we discussed above, the LHC signal of the model is multi-photon + missing or multi- 
lepton + missing. In the following, we list the search modes which we investigate and the 
set of kinematical cuts we used for the optimization of the significance. 
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We perform exclusive searches on the number of photons and of leptons. For the 
photon mode, we analyze one photon (I7) mode and two or more photon (2 > 7) mode. 
For the lepton mode, we analyze zero lepton (0£) mode , one lepton (1£) mode, same- 
sign two leptons (SS2£) mode, three leptons (3£) mode and four or more leptons (4 > £) 
mode. Each mode is divided into submodes according to the number of leptons (for the 
one photon mode) or the number of tau jets (for the lepton modes). These modes are 
summarized as: 

'l7 + (l£,2> 



Photon modes 



Lepton modes < 



2 > 7 

0£ + (2r, 3 > r) 
H + (Or, lr, 2 > r) 
SS2£ + (Or, lr, 2 > r) 
3£ + (Or, lr, 2 > r) 
4 > i 



(10) 



As mentioned in the previous subsection, we use a set of kinematical cuts for opti- 
mization of the significance. More concretely, we prepare steps of cuts for the missing 
energy, number of jets, jet p?, photon p^, lepton p^ and tau-jet px as follows: 



{^T.miss > 50, > 100, > 150, > 200, > 300}(GeV), 
{iVjets > 0, > 1, > 2, > 3, > 4} 

where p T (ji) > 100 GeV, prO'2,3,...) > 50 GeV, 
{ptU) > 0, > 100, > 200, > 300, > 500}(GeV), 
{p T (7) > 30, > 60, > 90}(GeV), 
{p T {£) > 10, > 20, > 30, > 50, > 70}(GeV), 
{p T (r) > 20, > 30, > 40, > 60, > 80}(GeV). 



(11) 

(12) 
(13) 
(14) 

(15) 
(16) 



We show in Table |2] which set of cuts are used for optimization in each mode. For 
example, in the three leptons + two or more tau-jets mode, the cut on the missing energy 
-^T.miss = |PT,miss| is selected from the set in Eq. ( ITT]) , cuts on the first two jets, ji and 
J2, are selected from Eq. ( TLBl . cuts on the first and third leptons, £\ and £3, are selected 
from Eq. f lT5|) and cuts on the first and second tau-jets, T\ and r 2 , are selected from Eq. 
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Tab 



e 2: Kinematical cuts used for each mode. 
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3.3 LHC Reach 

In Figs. [2] and [3j we show the discovery region of the gauge mediation model. 

Among the searching modes we analyzed, the most important modes are the two 
photon mode, same-sign two lepton mode, three lepton mode, and the four or more 
lepton mode. The results for these modes are shown in the figures. The red region is the 
region already excluded by the Tevatron trilepton [191 12D] and photon signal search [21] . 
The blue region is the region excluded by the LEP experiment [22] . 

In the figures, we show the discovery region for 1 fb _1 at the 7 TeV run and 1 and 
3 fb _1 at the 14 TeV run. First of all, we see a clear separation between the discovery 
regions of 2 photon mode and multi-lepton modes. As discussed above, this corresponds 
to a change of the NLSP shown in Fig. [TJ The result for the 2 photon + missing search 
is illustrated in the top left figures labeled with (a). The discovery region is determined 
by the gaugino masses, especially the wino mass. 

In the stau-NLSP region, the multi-lepton modes are effective. The three displayed 
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Figure 2: LHC discovery region for tan /3 = 10. The red region is the region already 
excluded by the Tevatron trilepton [191 120] and photon signal search [21] . The blue region 
is the region excluded by the LEP experiment 
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Figure 3: Same as Fig. |2} but for tan (3 = 40. 
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Figure 4: Combined discovery region tan ft — 10 and 40. The red and blue regions are 
same as in Figs. [2] and El The purple lines show the points realized in the minimal GMSB 
with the number of messengers N 5 = 1 to 5. On each line, the bold line shows the region 
with A < 80 TeV (see Sec. @J). 

modes, same-sign two lepton, three lepton, and four or more lepton modes have similar 
discovery regions, which are roughly determined by the mass of squarks, whose production 
is the dominant SUSY production. We should mention that other lepton modes, such as 
zero-lepton multi-tau mode and one lepton + multi-tau mode have narrower discovery 
regions than the above multi-lepton modes, even in the tan ft = 40 case. This is because 
a small fraction of produced tau leptons are detected as tau jets because of its low identi- 
fication efficiency and there are comparable amount of e//i leptons which come from the 
decay of the tau leptons. Thus even for high tan ft cases, the multi-lepton modes are more 
important than the less-lepton + multi-tau modes. 

Finally in Fig. HJ we show the combined result for the discovery region. Discovery 
regions for all modes listed in Sec. 13.2.31 are combined. In the figures, we also illustrate 
with purple lines, which correspond to the minimal gauge mediation models with the 
number of messengers N$ = 1 to 5. 
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3.4 Comments on More General GMSB Models 



We have assumed a "GUT relation" for the gaugino and scalar masses. However, in 
general GMSB models, this relation is not always maintained. Although our present 
analysis cannot be directly applied to such models, the story does not change significantly, 
as long as the slepton or bino-like neutralino is the NLSP. This is because although details 
of the SUSY signal depend on each SUSY spectrum, the SUSY signal with multi-leptons 
+ missing energy or multi-photon + missing energy is naturally expected in the case of 
the slepton or neutralino NLSP for any SUSY mass spectrum. Roughly speaking, if the 
number of SUSY events is O(10 — 100), SUSY can be discovered, as shown in Fig. [51 In 
this figure, we show the scatter plot of the total SUSY production cross section a and 
the significance Zb for the integrated luminosity 1 fb _1 , assuming the "GUT relation" 
for the gaugino and scalar masses. One can see that in the bino NLSP case, almost 
all the region where the event number exceeds about 20-40 can be discovered. This is 
because, in the bino NLSP case, the discovery relies almost only on diphoton + missing 
energy and other objects such as jets or leptons are irrelevant. In contrast, in the stau 
NLSP case, the required number varies widely, O(10 — 100). In some parameter points, 
the SUSY cascade decays tend to emit tau-lepton instead of e//x leptons, or emissions 
of high-energy leptons are suppressed because of kinematical reasons. In such cases, the 
required number of SUSY productions becomes large. In contrast, if the SUSY spectrum 
prefers high-energy e/fi lepton emissions, the required number can be less than 50. It 
is expected that this argument can be applied to more general GMSB models except for 
specially-tuned parameter points, as long as the slepton or neutralino is the NLSP. Thus 
we expect the SUSY particles can be discovered if the number of the SUSY particles reach 
0(10-100). 

In Fig. [61 we show the masses of SUSY particles for which 100 pairs of SUSY particles 
can be produced for a given integrated luminosity. If the required number of SUSY events 
for discovery is 100, as discussed above, the lines in Fig. [6] corresponds to the reachable 
mass range. 
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Figure 5: Relation between the significance Zb and cross section a for the integrated 
luminosity 1 fb~ . Left: the stau NLSP case, right: the bino NLSP case. 




Figure 6: Reachable SUSY masses at the LHC as a function of the integrated luminosity. 
Assuming that 100 SUSY events are required for discovery, the lines show reachable SUSY 
masses for a given integrated luminosity. The left-hand side is for the 7 TeV run and the 
right-hand side is for the 14 TeV run. 
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4 Implications for Minimal GMSB Models 



In the previous sections, we have analyzed the LHC signatures and discovery potential of 
the LHC in a setup more general than the minimal GMSB models. In this section, we 
discuss implications of our analysis for the minimal GMSB models. 

Although, as described in the introduction, low-scale GMSB models are very attrac- 
tive, they tend to suffer from a serious problem: instability of the SUSY-breaking vacuum. 
This problem stems from the fact that once one introduces the messenger particle in the 
theory for mediating the SUSY-breaking effect, a supersymmetric vacuum, in which the 
messengers develop the VEV, often comes into the theory and the former SUSY-breaking 
vacuum then becomes only a metastable vacuum. Two possibilities can be considered to 
avoid this unwanted problem. One is to force the metastable vacuum to have a lifetime 
much longer than the age of the universe, and the other is to construct a model with a 
stable SUSY-breaking vacuum. 

However, the mass scale of the SUSY particles is strongly constrained in both cases. It 
is known that if the SUSY-breaking vacuum is stable, the gaugino masses are suppressed 
[23j 121]. In GMSB models with a perturbatively calculable stable vacuum, there is very 
strong upper bounds on the gaugino masses [25]. Combining the updated null- result of 
diphoton+missing energy search at the Tevatron [21], such types of GMSB models are 
completely excluded if m 3 / 2 < 16 eV. The remaining possibilities are GMSB models with 
a metastable vacuum or strongly interacting GMSB models with a stable vacuum [2"6"| [27] . 
In the following discussion, we concentrate on the case with a metastable vacuum. 

There is also a strong constraint for the MSSM mass spectrum in this case [5]. To 
achieve a SUSY-breaking vacuum, whose lifetime is longer than the age of the universe, 
there is an upper bound on A g and A s . For example, if the number of messengers N5 = 1 
and if we adopt the IYIT model [2H] with an SP(1) gauge theory as the SUSY-breaking 
sector, then A 9 must satisfy A g < 80 TeV [8] . 

This upper bound on A 9 (and therefore A s ) imposes severe constrains on the MSSM 
mass spectrum. For example, if the number of messengers is one and tan/3 is not so large, 
then the NLSP is the neutralino and the region A g < 124 TeV has been already excluded 
by the Tevatron diphoton+missing search [21] as shown in Fig. H] (a). 
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Figure 7: Required luminosity for testing the minimal GMSB with N 5 = 1,3, and 5 for 
varying tan /3. 



One possible way to evade the above experimental bound is to increase the number of 
messengers N$. However, there is also an upper bound on the number of messengers N$ 
to realize a successful GUT unification, which is an important motivation for introducing 
low-energy supersymmetry. A straightforward one-loop calculation yields the condition 

< 5. According to a more detailed analysis in Ref. [29], the necessary condition further 
strengthens to N 5 < 4 for M = (9(100) TeV, otherwise the predictable GUT unification 
is spoiled. Combining this constraint with the one resulting from the metastability con- 
dition, we have constraints on the scales A g < 320 TeV and A s < 230 TeV. 

In Fig. HI the purple lines indicate the minimal GMSB models with Ns = 1 to 5. On 
each line for N$, the bold region corresponds to the region with A g /N^ < 80 TeV, where 
the condition from the stability of the SUSY-breaking vacuum is satisfied. We see that 
for N5 < 4, the allowed region can be tested at an integrated luminosity 3 fb -1 at 14 TeV 
for both cases with tan /3 = 10 and 40. 

In Fig. [7J we show the integrated luminosity required to test the minimal GMSB 
models with A g /N§ = 80 TeV for varying tan/3 and N$ = 1,3, and 5. The NLSP is 
the stau except for the region tan/3 < 35 with A^ 5 = 1, which is already excluded by the 
Tevatron experiment. Therefore, the relevant discovery mode is multi-lepton + missing 
energy. We note that when N 5 = 1, almost all the regions in tan/3 are already excluded 
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by the current experimental limits, but a small region above tan /3 > 35 survives. This 
small region can be tested at a very early stage with 1 ftT 1 at 7 TeV. We also see that 
even for N 5 = 5, all the regions in tan/3 can be tested with an integrated luminosity 
of about 10 fb _1 at 14 TeV. For larger tan/3, larger integrated luminosity is required to 
test the models. The reason is that more tau leptons instead of e//x leptons tend to be 
emitted in SUSY cascade decays for larger tan/3, because of both kinematics and the 
SU(2)l interaction carried by the mixed stau. 

5 Conclusion and Discussion 

In this paper we have investigated the LHC discovery region for low-scale GMSB models 
with a setup applicable for a wider class of models other than the minimal GMSB models. 
We have performed a comprehensive study while giving careful treatment for realistic 
experimental setups. Our result is thus a very conservative one. Although we have set 
the lepton and photon detection efficiencies are 100 % for the SM background, if we 
impose the same fake rates as used for the SUSY signals, the reduction of the background 
is expected, which leads to a wider discovery region. In addition, we have used the leading- 
order SUSY production cross section. Generally, the next-leading order cross section is 
larger than the leading-order one. Thus, if we adopt the next-leading cross section, the 
discovery region is extended by about 10 %. 

We have applied our result to the minimal GMSB models, in which the MSSM SUSY- 
particle masses have strong upper bounds from the stability condition of the SUSY- 
breaking vacuum. We have shown that all the region can be tested at an early stage of the 
LHC if m 3 / 2 < 16 eV. Let us comment on the case with a lighter gravitino. Approximately, 
the upper bound of the scales A g and A s is proportional to y/m 3 /2. If in the future the 
upper bound on the gravitino mass will be improved to, say, 3 eV, then it will be predicted 
that the SUSY can be discovered at a very early stage ^/s = 7 TeV, even if N 5 = 5. 
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